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Abstract 

Gluon distributions in real and virtual photons are calculated using evo- 
lution equations in the NLO approximation. The quark distributions in the 
photon determined on the basis of the QCD sum rule approach in ref.[l] are 
taken as an input. It is shown that gluon distribution in the photon can be 
reliably determined up to x = 0.03 -j- 0.05 much lower than the corresponding 
values in the case of quark distributions. Two variants of the calculations are 
considered: (1) it is assumed that there are no intrinsic gluons in the photon 
at some low normalization point Q 2 = Ql ~ lGeV 2 ; (2) it is assumed that 
gluonic content of the photon at low Ql is described by gluonic content of 
vector mesons p,cu,cp. The gluon distributions in these two variants appear 
to be different. This fact permits one to clarify the origin of nonperturbative 
gluonic content of the photon by comparing the results with experiment. 

Structure functions ^(x) for real and virual photon are calculated and it 
is shown that in the region x > 0.2 where QCD approach is valid, there is a 
good agreement with experiment. 
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1. Introduction 



In paper [1] the structure functions of real and virtual photons were cal- 
culated. The hadronic part of the photon structure function was calculated 
in a modelless way on the basis of the QCD sum rules applying to structure 
functions [2,3]. In this aspect the calculation made in [1] differs from earlier 
considerations of the same problem [4-7]. The idea of the ref.[l] consider- 
ation was the following. First, the structure function of the virtual photon 
(photon-target) was studied under conditions when photon virtuality p 2 < 
and | p 2 |>> R~ 2 , R c - is the confiniment radius, but at the same time 
| p 2 \« Q 2 = — q 2 , where q - is the photon projectile momentum. In ex- 
pansion in power in p 2 /q 2 only terms of the first nonvanishing order were 
taken into account - the terms of the lowest twist-2. The operator product 
expansion (OPE) over 1/p 2 was performed and the first and second order 
terms were taken into account. The contribution of the latter corresponds 
to interaction of quarks with vacuum gluonic fields and leads to appear- 
ance of terms proportional to gluonic condensate in the structure function 
of the virtual photon. Thus, the structure function of the virtual photon 
F2(s,p 2 ), s = (p + q) 2 was constructed as a series in 1/p 2 . On the other 
hand, using analytical properties of F2(s,p 2 ) as functions of p 2 by means of 
dispersion relations in p 2 at fixed s F2(s,p 2 ) was presented through contribu- 
tions of physical states modelled as contributions of vector mesons (p, ip, u) 
and continuum. By identifying both of these representations the structure 
function of the virtual photon including its hadronic part was found. At such 
approach the structure function possesses correct analytical properties in p 2 
and has no fictitious singularties at p 2 = inherent to bare diagram of Fig.l 
(for massless quarks). 

The result of the calcultions of the transverse photon structure function 
is [1]: 



(— Fj(x,p 2 , « ! )^Ee!{[2 + «(x)(/n V 2 , - 3+ 

7T q H X Z {So—p Z ) 

So— V 2 {p — rriy q ) z 27sqX z tt p 

where x = Q 2 /2v is the Bjorken variable, k(x) = x 2 + (1 — x) 2 , sq = l.bGeV 2 
is the continuum threshold, my iq are the vector meson masses, p and u for 
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q = u,d and (p for q = s. (We restrict ourselves to consideration of three 
flavours) . 

As has been shown in [1] and as is seen from (1) the first term in square 
brackets in the right-hand side of (1) corresponds to the continuum contribu- 
tion in dispersion representation in p 2 while the second to the vector-meson 
contribution. Applying duality considerations we arrive at the conclusion 
that the first term in the right-hand side of (1) can be considered as a pho- 
ton perturbative contribution into F^ and the second term as hadronic part 
of the structure function F<i- Thereby, we unambiquously select from the 
diagram of Fig.l the soft and collinear quark contribution which can referred 
to quark sea distribution in photon but not to direct contribution of photonic 
operators in OPE. 

At x > 0.2 the gluonic condensate contribution in (1) is less than 15% 
even at p 2 = 0. This allows one (at these values of x) to extrapolate (1) to 
the point p 2 = and to get structure function of the real photon. 

An essential disadvantage of the discussed method of the structure func- 
tions calculation [2,3,1] is that it is inapplicable at small x and x close to 1. 
As is seen from (1), the correction term of OPE proportional to gluonic con- 
densate fastly increases with x decreasing, i.e. the series in 1/p 2 diverges at 
small x. This is caused by the fact that in the imaginary part of the forward 
7 — 7 scattering amplitude - the diagrams of Fig.l - the quark virtuality at 
vertical quark lines k 2 ~ xp 2 for massless quarks. Thus, at small x quarks 
appear to be in near the mass shell, i.e., in the nonperturbative region. At x 
close to 1, the situation is similar [3]. As was shown in [1], eq.l is true for the 
real photon at 0.2 < x < 0.7. With | p 2 | increasing the applicability region 
of (1) expands. 

In expression (1) for the structure function F2(x,p 2 ,Q 2 ) the evolution of 
F2(x,p 2 , Q 2 ) with Q 2 was not accounted, because for the solution of the evo- 
lution equations it is necessary to know F2 throughout the whole interval 
of x. That it why one may expect that (1) is true at intermadiate values 
Q 2 ~ 5 — lOGeV 2 , where the role of evolution is still inessential. (When com- 
paring with the experiment it turned out that (1) well describes experiment 
even for Q 2 « 20GeV 2 ). 

In this paper we will try to partially correct this disadvantage, to take into 
account the evolution equations and to calculate gluonic distribution in the 
real and virtual photon. The calculation of gluonic distribution is preferable 
by virtue of our approach because gluons are produced by quarks and gluon 
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energies are essentially smaller than quark energies. This allows us to hope 
(and is confirmed by calculation) that basing on the known quark distribu- 
tions in the photon at x > 0.2 we shall be able to find gluonic distribution 
at noticeably smaller x. (The x > 0.8 region weakly affects gluonic distribu- 
tion at x < 0.5). The interest to gluonic distribution in the photon is also 
connected with the proposal to measure it on HERA accelerator [8]. 

Calculations will be performed in the leading (LO) and next to leading 
(NLO) approximation. We consider two variants of the boundary conditions. 

The first variant is based on an assumption that all gluonic distribution 
in the photon stems only from their emission by quarks and that the photon 
has no intrinsic gluons of nonperturbative origin. I.e., we will believe that 
at some low normalization point Q 2 = Qq the gluonic distribution (as well 
as the sea quark distribution - we determine them in Sec. 3) vanishes. This 
idea is close to the forwarded earlier analogous idea on the origin of gluons 
and sea quarks in hadrons (see, e.g. [9]). It seems that in the case of a 
photon, espesially of a virtual one, such an idea has much more rights for 
existence than in the case of hadrons. Let us once more emphasize that the 
main difference between this consideration and the previous ones (see, e.g. 
[10-20]) is that the hadronic part of the structure function of a photon at 
Q 2 = Qq is taken from [1] and the point-like photon is separated from its 
hadronic components in a way as it was explained above. 

The second variant of the boundary conditions is also related to the results 
of [1]. The second (hadronic) term in eq.(l) corresponds to the contribution 
of vector meson, i.e. to transition of a photon into vector meson with sub- 
sequent scattering of the projectile virtual photon on this vector meson (see 
Fig. 2). Therefore one may speak about gluonic sea of the photon which corre- 
sponds to the vector meson intrinsic gluons (naturally, with a corresponding 
normalization). 

As follows from our results, the gluonic distribution appears to be essen- 
tially different in these two variants. Thus, comparison of gluonic distribution 
with experiment would clarify hadronic structure of photon. 

In Sec. 2 we discuss the evolution equations in the NLO approximation. In 
Sec. 3 we will find their solution and discuss the results obtained. 
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2. Evolution Equations. Partial Solution 
of Inhomoqeneus Equation. 

As has been originally shown by Witten [10], evolution of the photon struc- 
ture function differs from the standard evolution of the structure functions 
of hadrons. The reason for this is that in the photon case there exists an 
additional mechanism of direct (point-like) production of quarks from pho- 
ton (as well as of gluons in the next order in a s ). This mechanism results 
in logarithmic increasing of the structure function with Q 2 unlike evolution 
of hadronic structure functions where quarks and gluons may be produced 
only from quarks or gluons. Formally, in the OPE technique the moment 
expansion of the photon structure function has the form [10] 



F?(Q') = J Q dx-x 



n-2 



F*(x, Q 2 ) = E C n (Q 2 /^ 2 ,g^ 2 )) < 7 | 0\ 

i=NS,S,G 



7 > + 



+ Cl{Q 2 /v\g{i?))< 1 \Ol\ 1 > 



Here 0\ = (Of ,0^, O n 



G\ — 



(2) 

O n , correspondingly, nonsinglet and singlet 
quark and gluon operators of spin n and twist-2 and O 7 are the photon 
operators (see e.g. [11]), C l n (Q 2 / fj 2 , g 2 (//)) is the column of Wilson coefficients 



C n (Q 2 /fi 2 ,g(fi)) = 



= C n (Q 2 /fi 2 ,g(n)) 



C^(Q 2 /fi 2 ,g(fi)) 
C s n {Q 2 /v 2 ,g^)) 
C^QV/i 2 ^/,)) 

As was noted in [10], the second term in (2) for the case of hadrons would be 
the 0(a) correction, since C^(Q 2 / fj 2 , g(fi) ~ 0(a). For the case of photons 
this term should be also taken into account because matrix element < 7 | 
O l n I 7 > in the first term of (2) is also of order a (< 7 | O^ | 7 > is, naturally 
of order 1). 

Thus, the renormalization group equation for C n takes the form [10]: 



C n 


( In 




C n 


C 7 

n 




0) 


C 7 

n 



(3) 



where % is the standard matrix of anomalous dimensions and K n = 
(K^ s , K%, K„) leads to mixing of photonic and hadronic operators. One 
can find the explicit form of % and K up to terms aJnQ 2 in refs. [12,13]. 
(See also expressions 8-10 ). The term K n gives an additional contribution in 
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the Gribov-Lipatov-Altarelli-Parisi (GLAP) evolution equation to the quark 
distribution function from the direct quark and gluon production by photon 
(see Fig. 3). Q Solving eq.(3) and substituting the result into (2) we can find 
the final expression for ^(Q 2 ). In symbolical notation it looks like (see, e.g. 
[10, 12]): 

F?(Q 2 ) = 4(1, f)M n {g\ 9l) < 7 I O n {Ql) | 7 > + (4) 

Mf,g 2 o)Cn(hg 2 ) + cz(i,g 2 ) 

where 

f = 47ra s (Q 2 ); g\ = 4 7 ra s (Q 2 ) ); 



So 

M n (g\gl) = exp J dg 
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in{g) 



X = ]d 9 J ^.M"(f,g>) 
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Functions 



°- (1 '^- (1+ 16^^ >'[ <e 2 > for S 

C^f)=<e*>^- 2 B G n (5) 

CKhf) = ^ 3 £ A Bl- < e 2 >= (1//) £ e 2 
4tt i=1 i=1 

£7 = (2//) and 5^ are presented in [12]. 

In the general case X n , M n , A f n ,~K n - are matrices. In what follows all 
calculations are made on an example of nonsinglet part of : F^ s . In this 
case M n , j n , K n are algebraic expressions. (The calculation of the singlet part 
is principally the same.) 

For the nonsinglet case (4) is reduces to: 

F n ( NS) = c xs {1>f> {M " NS (f,9$) < 7 I CWo) I 7 > +X» s (?,gt)} (6) 



x It should be reminded that K play the role of anomalous dimensions and are p 2 independent. All p 2 
dependence in (2) is gathered into matrix elements < 7 | O l n | 7 > and appears in fact as a boundary 
condition. 
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The factor in the curly brackets is usually treated as the moment of distri- 
bution function (see, however [21], when the other difinition of the distri- 
bution functions was discussed, in which C^ s also was included in it). In 
the leading logarithmic approximation C l n = 1 and the standard definition 
F% = Eefg^, q l n = J x n ^ 1 qi(x)dx is reproduced. In NLO the moments of 
quark distribution function q NS are defined by 



«r = / X- 1 q NS (x)dx = M« NS (?,g>) < 7 I N n s | 7 > +X» s (g\ gt) 

In the same way the moments of singlet and gluon distribution functions 
are defined. The first term in (4) corresponds to the contribution from 

S,S,G (Q 2 / /J? i g 2 ) in (2) (with account of the evolution (3)), the second 
and the third terms arise from C^(Q 2 /ii 2 ,g) in (2) as a result of solution of 
eq.(3). Usually, the first term in (4) is treated as hadronic part of the struc- 
ture function, since its form coincides with the solution of renormalization 
group equation for hadrons; the second term - as "point-like" production of 
quarks and gluons by photon (which involves characteristic In Q 2 behaviour) 
and the third term - is the proper photonic part which is included into quark 
distributions in some factorization schemes (e.g., DIS~ P in ref.[16]). How- 
ever, such separation may be dangereous because of the possibility of double 
counting, as was discussed in [1]. One can make this separation correctly 
basing on approach and results [1] where it is performed in the framework 
of QCD sum rules. To this end it is more convenient to use the evolution 
equation language. 

The explicit form of the evolution equations for the distribution functions 
qi(x) and G{x) was given in [14] in the leading logarithmic approximation 
while the next order terms were taken into account in [11,15] . These equations 
have the form: (i means flavours of quarks) 



dG ° (P gq e £(<f + q l ) + P gg e G) + k g e r 



dt " 27T ' H > ' *99^"J ' 2yr 

where t = In Q 2 /A 2 ; P = P° + as< ^ ^- P 1 , and P-are the splitting functions 
[11-13], 

K* = {Kf + ^-Kf) K G (x) = ^ Kh(x) (7a) 
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Here P q stands for convolution 



f —P(x/y)q(y,t), 
x V 

T(y) = 8(1 — y) is the photon distribution function in the photon. 
a/2nK qi G{x) are related with anomalous dimensions K n = (K^ , K^) 
by Mellin transformation. The exact relation can be found in paper [11]. We 
do not present it here since in what follows we shall need only K q ^(x), the 
explicit form of which is [11] : 



Kf' l \x) = 3e 2 ^(x), i = u,d,s 



(8) 



where k®(x) = k(x) is defined in (1); 



K i(x) = -e-{4 - 9x - (1 - 4x)lnx - (1 - 2x)ln 2 x + 4/n(l - x)+ (9) 

+ [4/nx - 4lnxln(l -x) + 2ln 2 x - 4/n(l -x) + 2/n 2 (l -x)- (2/3)tt 2 + 1Q]k(x)} 

Kq(x) was taken from [16] where the error in [12] and in the papers based 
on it was corrected 



K h( x ) = Ee?-4[-16 + 8x + ^x 2 + ^- (6 + 10x)/nrc - 2(1 + rc)/n 2 rc] (10) 

Here, analogously to (6), the function F2 is determined through quark and 
gluonic functions as [15, 11] 



F 2 (x) = x\J2 e{(gi(x) + 4^ ® «)+ 



+ < e 2 > —Eg G + 3 E ej5 7 (x)a/47r 

47T 



where q = q + q 



(ii) 



= (8/3) { ^ - - \ 



1 + x 2 
1 — X 



+ 



+ 



+ (l + x 2 ) 



/n(l — x 
1 — x 



- (9/2 + tt 2 /3)(5(1-x) 



J + 



B G (x) = 2f - 2x + 2x 2 ) In— 1 + 8a: (1 - x) j ; 

B,(x) = (2/f) ■ B G {x)) 

B% and B^ denned in (5) are n- moments of these functions. (All expressions 
here are in the MS scheme). The last terms in eqs.(7) correspond to the 
process of point-like (direct) production of quarks and gluons from photon 
in LO and NLO: we preserve only the terms proportional to InQ 2 or a s lnQ 2 
(see Fig. 3). 

In the case of DIS 1 scheme suggested in Ref.[16] 

Ff'^x) = x k e? (f'*r(z) + £ Rt e q DIS ->(x)) + 

+ < e 2 > ^B G (SG(x)\ (11a) 
4tt J 

where the last term in (11) is included in the q DISl (x). This redefinition leads 
to transformation of K^ G (= K 1 ^ 1 or K G , see 8-10). 

K (i) dis, _ K {i) 

where SK qG from eq.(3.1, A.11-A.12) in ref.[16]. This scheme has many 
advantages in comparison with the MS scheme (see discussion in [16]). 

Especially in the DIS 1 scheme there is no strong difference between LO 
and NLO results. Hereafter we will work with the DIS 1 scheme (and omit 
DISj indeces). We would like to stress, however, that our results are com- 
paratively insensitive to the choice of renormalization scheme (MS or DIS 7 ), 
see discussion at the end of the paper. 

The evolution equations (7) are equivalent to the renormalization group 
equations for the moments (3). The form of equations (3) and (7) is inde- 
pendent of the photon virtuality, the p 2 dependence is hidden in the matrix 
elements < 7 | O n \ 7 > and therefore appears only through the boundary 
conditions. Evidently, for virtual photon the renormalization group equa- 
tions are valid at Q 2 >> — p 2 . In what follows we consider only the case of 
transversally polarized real or virtual photon and for simplicity the notation 
F2 will be used for at p 2 7^ 0. The unpolarized virtual photon case will 
be discussed at the end of Section 3. 

The evolution equations (7) are inhomogeneous. As is known, the general 
solution of inhomogeneous equation is given by the sum of partial solution of 
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inhomogeneous equation plus general solution of homogeneous equation. The 
latter coincides with the standard solution of evolution equation for hadrons. 
We choose the partial solution of evolution equation in NLO by requiring 
its vanishing at Q 2 = Qq (see refs. [15-17]). The solution of inhomogeneous 
equations, represented in terms of moments of singlet and gluon distributions 
is given by [16] 





a 2 f 




~~ 2^ Jo I 



A 



+>—\_p_ + 

1 _ £l-2A+//3 



2n/a s P+ - 

P-AA 



A) 



+ 



1 - 2A+/A, 
2 



A 



+ -X-- A)/2 

1 _ ^-2A+//3 



X 



2A+/A) 

A) 



P+ #1 " W^P+Ko- 



A 



P+P\P+ - ^ A + P + - 



P+Pip 



2 A_-A+-A,/2, 



+ (terms with interchange A + — > A_, A- — > P 



)} 



+ 



Here A, A-? P_ are matrices; = E;=i (? n + # 



in 



(12) 



P = 



pi pi 

99' 9.9 

pi pi 

59 55 



P + = ±- 



Pn 



A+-A_ 



Pn = 



pO pO 

99 9.9 

pO pO 

.99 53 



A ± = l/2(i* +i* ±^(i* 



33 



pO \2 _|_ 4p0 pO ^ 

33^ 1 93 39' 



99 



P° = 2fP' 

93 J < 



93' 



#1 = 




2/ 



^(^ 2 ) . if (o,i) = v K l{ W 

<Qir q ^ q 



/ \ IV GJ u sWoJ 1=1 

Pqq 1 , Pqg 1 etc. mean the n-th moment of the standard splitting functions. 

As mentioned in refs. [15,16], such a choice of the partial solution when 
even at large Q 2 the terms L~ 2A ±/^ are retained permits one to get rid of 
nonphysical poles at n = 1.596 and n = 2 corresponding to 1 — 2A+/A) = 
and A + = 0. These poles would correspond to the poles in x-distribution of 
the type 1/x 1596 and 1/x 2 . The vanishing at the same points numerators 
result in the logarithmic divergences in x-distributions instead of power-like 
ones. The partial solution chosen in this way is completely analogous to 
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the term X(g 2 ,gl) in the solution of eq.3. As already has been mentioned, 
this term incorporates the terms proportional to InQ 2 , related to the point- 
like mechanism of quark and gluon production by photon in LO and NLO. 
However, it would be erroneous to think that the solution of inhomogeneous 
equation with zero boundary conditions takes into account the whole con- 
tribution from the process of the point-like quark production at a given Q 2 
while the general solution of homogeneous equation (equivalent of M n in (4)) 
- only hadronic part, i.e. the contribution of soft and collinear 
quarks Q. At such a choice, there may appear a double counting (or a deficit 
if we take hadronic part to be zero [16]). 

One of the methods to solve this problem correctly was discussed in the 
Introduction. This is to use as boundary conditions the results of ref. [1] - 
eq.(l),, where the point-like and hadronic contributions into photon structure 
functions - respectively, the first and the second terms in curly brackets 
in (1) - are calculated in QCD. In such case the mechanism of the point- 
like production corresponds to solution of inhomogeneous equation with the 
boundary conditions from the first term of (1) to which one should add the 
hadronic part corresponding to solution of homogeneous equation with the 
boundary conditions from the second term of (1). 

Owing to linearity of the evolution equations, it is, naturally, more conve- 
nient to take the solution of inhomogeneous equation with the zero boundary 
conditions and to add the solution of homogeneous equation with the bound- 
ary conditions which are the sum of the boundary conditions of the point-like 
and hadronic parts. But in doing so, the solution of homogeneous equation 
is not only the contribution of hadronic part (e.g., it is, erroneous to choose 
in this case the boundary condition from VDM considerations). 

Before going to solution of homogeneous equation, let us come back once 
more to the general form of inhomogeneous equation (7). 

Consider two functions cp(x) — q NS (x) = (l/2)(u(x) — d(x)) and r = 
—(l/2)(u/e 2 l — d/e 2 i ) where is the charge square of the corresponding quark. 
The evolution equation (7) and the boundary condition on are the same 
for d and s quarks besides small difference in the boundary condition (1) 
arising from the mass difference of p and (p mesons. (The calculation shows 
that the influence of this difference on gluon distribution is very small) 

Thus, we may safely take d(x) = s(x), and, correspondingly, 

2 Such a distinction may be achieved, as we shall see in what follows, but only at certain boundary 
conditions. 
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d(x) = s(x) : ft = + =4^| + ^ (13) 

where A = e 2 u — e 2 d = 1/3. The inhomogeneous evolution equation (7) ex- 
pressed through these functions is: 

^-^(eJ-4)3k(.) + ^)(.)) + 



ttsW '*-P„(*/»)rtv,t) 



2i ; s 

dG{x,t) _ aa s (t) x a s (t) j 



4Ee 2 

^G(«/j/)(-^ L ¥'(y»*)+ 



4/f ''^-r(2/,t))+P 5 ,( a :/|/)G( 2 /,t) 



A 



V 



where notations are the same as in (7), (7a), (8). The evolution equation 
for t(x) has the form of standard evolutions for hadrons and all quark and 
gluon point-like production is now absorbed in <p function. Thus, <p(x) may 
be treated as valence quark distribution in the photon and A(e 2 u e 2 d / £s)t(x) as 
"sea quark" distribution f\. In the next Section we consider two variants of 
boundary conditions and discuss the obtained solutions. 

3. Solution of Evolution Equations. Discussion of Results. 

In the first variant of the boundary conditions we assume that at some, 
relatively small Q 2 = Ql there are no r ("sea") quarks, as well as gluons but 
there are only "valence" quarks ((^-quarks). 

Physically, this corresponds to assumption that both gluons and sea quarks 
in photon are produced in perturbative way due to their emission by valence 
quark (as well by a bare photon for the gluon case) and they are absent 
at the low normalization point Q 2 = Ql ~ lGeV 2 . As has been already 
noticed in the Introduction, an analogous idea was forwarded to describe 
the structure function of hadrons where it was supposed [9], that at small 
Q 2 = Qo ~ IGeV 2 a hadron consists only of valence quarks. But such an 

3 Such detemination of "valence" and "sea" quarks distribution is close to that used in [20]. 
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idea failed in description of the hadron structure functions - it appears to be 
inconsistent with experiment. This idea has better chances for photon since 
it can be naturally thought that photon has no gluons and sea quarks (in 
the sense defined at the end of the previous Section) as their constituents. 
Note that for virtual photon Q 2 should be chosen larger than | p 2 |. As 
the hadronic component contribution into structure function (1) decreases 
rapidly with | p 2 | increase, and hence it follows that for strongly virtual 
photon our requirement is fulfilled automatically. 

The total solution for the quark distribution function consists of the sum 
of solution of inhomogeneous equation which was under discussion in the 
previous section and of the solution of homogeneous equation 

<MM)_a.w } [Pm{x/y) ^ t) /y 



dt 2tt J ' r > ' " y 

dr^tl = a^)j tM + l Pw(a;/2/) _|_ G (y,t)& (15) 

with the boundary conditions from (1). Since we restrict ourselves to NLO 
of inhomogeneous equation i.e., to terms a s lnQ 2 , then in the same order the 
homogeneous equation can be solved in the leading logarithmic approxima- 
tion, so = j-j.. In this approximation in eq.(lla) for ^(x, Ql) (note, 
that G{x, Ql) = 0°) 

(1/x) F 2 (x, Ql) = E e? [ qi{x) + ^ qi (x) } (16) 

one should omit the term (a s /47r) © qi , since the solution of homogeneous 
equation qi(x,Q 2 ) has no InQ 2 and the contribution from the solution of 
inhomogeneus equation at Q 2 = Ql is zero. As the estimates show, this 
term gives a small contribution to the gluonic structure function at x < 0.5. 
Therefore the boundary conditions have the form 

F 2 (x,Ql)=xj:e 2 qi(x,Q 2 ) (17) 
G(x,Ql) = r(x,Ql) = 
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Accounting for (1), (13) we have 



t ) = ( [2 + 9 ° - 3 + -^)]+ (18) 

47T I X 2 (Sq — J9 Z ) Sq — V 



4 pM 



2{rriy — p 2 ) 2 
where to = ln(Ql/A 2 ) 

and my means the weighted average of p and mass. Let us represent ip(x, t) 
as a series 

*) = t ) + £ ^M(/nf ) m (19) 
i m! t 

For r(x,t) and G(x,t) similar expressions take place: 

r(x,t) = t Tj ^-(lnU m , G(x,t) = t^-(ln-T (20) 
i ml to i m! to 

Substituting (19), (20) into the evolution equations we get recurrent relations 
(/ = 3) 

fcW = ( J-r K, ® < ® --^ t ) 

Po x . ' 

m times 

n(x) = o 



2 

C m+ i(x) = — 
Po 



X 51 e lPgq <Pm(%) H T" ~~ Tm(^) © ^ + © G m {x) 



2 r „ n . , , 1 A 



Po e d 

In order to solve equations (21) it is necessary to know the function ip(x, Qq) = 
(p(x,to) in the region 0.75 < x < 1 where eqs.(l) and (18) are invalid. Notice 
that the uncertainty at large x arises only in the second term in (1) or (18) 
related to the vector meson contribution. The natural extrapolation comes 
from quark counting rules and corresponds to f p (x) ~ (1 — x) at large x. We 
considered also the other forms of extrapolation f p in the domain of large 
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x: f p {x) = A(l — x) a , where A was determined from the matching with the 
second term in (1) at the point x = 0.75. It was found that G(x,Q 2 ) at 
x < 0.5 varies only by a few per cent if a varied from 1 to 2. We emphasize, 
that the account in the evolution equations of singular terms, proportional 
to 5(1 — x), 5'(1 — x) etc, (see [1]) would be a mistake. These terms (e.g. 
the term, proportional to operator a s ip I\ ip ip Ti ip) formally appear in the 
OPE over the soft photon virtuality 1/p 2 . However, as was explained in ref.3, 
these terms are unphysical because they correspond to the bump at x = 1 
and must be compensated in the correct theory. 

The series (19), (20) rapidly converge at x < 0.6 and Q 2 < lOOGeV 2 , so 
it is sufficient to take into account three terms in the expansion. In the 
homogeneous equation we can put the one- loop expression for a s : a s = 
4-7T / Poln(Q 2 1 A 2 ) , A = 230MeV, since this equation is treated in the lead- 
ing logarithmic approximation. In the case of inhomogeneous equation (10) 
the calculations are performed in NLO in DIS 1 regularization scheme and, 
correspondingly, two-loop expression for a s is used: 

J_ (n2] = 1 gi In ln(Q 2 /A 2 ) 

47i asW ] (3 Q ln(Q 2 /K 2 ) P'i ln 2 {Q 2 /A 2 ) 

The functions q(x), G(x) (q(x) = q s (x) is the singlet quark distribution) were 
obtained from the solution of the inhomogeneous equation (12) by the inverse 
Mellin transformation. The results are presented in Fig. 4-8. Dotted lines 
correspond to gluon condensate contribution, thin solid lines - to hadronic 
part and thick lines - total gluon distribution x G(x)/a. For p 2 ^ (Fig. 5-^7) 
thick solid lines correspond to transversally polarized photon case, and dashed 
lines - unpolarized photon case, which will be discussed at the end of the 
paper. 

As is seen from Figs. 4-7 the negative gluonic condensate contribution is 
essential only at x < 0.1, it increases steeply with x decreasing and for the real 
photon at Q 2 = 5GeV 2 comprises a half of the total G(x) at x « 0.02. This 
means that the QCD sum rule approach prediction for the gluon distribution 
function in the real photon is reliable up to x > 0.02 — 0.05 (in the framework 
of the chosen boundary condition - variant I (HI)). The gluonic condensate 
contribution decreases rapidly with increasing of photon virtuality —p 2 . 

As was discussed in the Introduction, such an essential extension of the 
applicability region towards small x is possible for the gluon distribution 
function only. For quark distribution the applicability region starts at larger 
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x > 0.1. 

As is seen from Figs. 4-7, the hadronic part is rather large for the real 
photon (~ 50% at Q 2 = 5GeV 2 and 30% at Q 2 = 50GeV 2 ) and decreases 
rapidly with | p 2 | increasing. 

An approach similar to ours for the case of the virtual photon was devel- 
oped in refs. [15,17]. Here, however, the hadronic part was not singled out 
from the boundary conditions and considered separately, as we did. 

The dashed lines in Fig. 4 at Q 2 = 5 and 50GeV 2 represent the results of 
ref.[16], where the NLO evolution equations were solved with the zero bound- 
ary conditions for inhomogeneous as well as for homogeneous equations. As 
is expected, they are much more less than our results although relative differ- 
ence decreases with x increasing. The difference is large even at Q 2 = 50GeV 2 
what means that the domain of Q 2 ~ 50GeV 2 is far from asymptotics, where 
dominant are the terms in the solution of inhomogeneous equations (~ InQ 2 
and a s lnQ 2 ) independent of the boundary conditions. 

Our results for the gluon distribution function are weakly dependent on 
the choice of the normalization point Qq. For the real photon this can be seen 
from Fig. 8. With Q 2 or | p 2 | increasing the Q 2 dependence becomes even 
weaker. The LO results are not too much different from NLO: the difference 
is less than 10-15% almost everywhere except for the highest considered vir- 
tually | p 2 |= 2GeV 2 where it is about 30%-40%. For this reason the LO 
results are not shown in the Figures. 

We turn now to the variant II (H2) of the boundary conditions discussed 
in the Introduction. The second term in the photon structure function, eq.l, 
corresponds to the target photon into vector meson transition with subse- 
quent projectile virtual photon- vector meson scattering. At such an interpre- 
tation it is natural to assume that at low Q 2 there are nonperturbative gluons 
in the photon and their distribution is determined by gluon distributions in 
vector mesons. As was shown in ref.[l], the second term in (1) is related to 
the transverse structure function ff (x) of the fictitious vector meson V built 
from one flavour quark q with unit charge 

(22) 

where p(x) is given by (18) and g\ is the 7 — V transition coupling constant. 
(Experimentally, gy/^n = 1.27.) fy(x)/x is equal to quark distribution, 
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fy(x) = x[q(x) + q(x) 

In a similar way we can define the gluon distribution Gy(x) in F-meson. 
Assuming the standard form of x-dependence, xGy{x) ~ (1 — x) ? \ we have 

x Gy(x) = c A 4 9 2 V {1 ~ xf (23) 
8tt z my 

where the normalization constant c is determined from the requirement that 
gluons in ^-meson carry 40% of y-meson momentum, c = 0.42. The gluon 
distribution in the virtual photon is related to xGy(x) by 



x G^(x,Ql) — — 2~E e 2 2 v x Gy{x) - 
9v \ m v ~ V ) 



^^-^-^wi-*) 3 ( 24 ) 

27r q [m v — p z Y 

Eq.24 is used as the boundary condition for the gluon distribution in the 
variant II of boundary conditions. In the same way the boundary conditions 
for sea quark distributions can be also imposed. But since they carry only 
10% of the vector meson momentum, their contribution to the final gluon 
distribution is negligibly small. We emphasize that (24) is the boundary 
condition for the hadronic part of the gluon distribution function, for non- 
hadronic, point-like part the boundary condition is still zero. As before, the 
boundary conditions (17) are used for valence quark distributions. (The ac- 
count of the order a s gluonic contribution to Fi{x, Qq) in (17) is beyond the 
accuracy of our calculation and, as the estimate shows, this contribution is 
small numerically). 

The results of the calculations in the variant II are presented in Figs. 9-12. 
All notations are the same, as in Fig. 4 8. In this variant the hadronic 
part dominates in the gluonic contribution function for the real or virtual 
photon even at Q 2 = 50GeV 2 . In comparison with the variant I G(x,Q 2 ) is 
essentially larger, especially in the case of the real photon - by a factor of 
2-3. The relative difference of G(x, Q 2 ) in variants I and II decreases with Q 2 
increased as a result of increasing role of point-like gluon production. The 
difference in gluon distributions in variants I and II is going down steeply 
with increasing of the photon virtuality | p 2 |. 

The results of ref.[18] at Q 2 = lOGeV 2 are shown in Fig. 13 by dotted 
line. The method of the calculations of ref.[18] is different from ours in three 
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aspects: i)low normalization point in inhomogeneous differential equation 
with the zero boundary condition at Q 2 = O.SGeV 2 in [18] comparing with 
our Qq = lGeV 2 ; ii) our boundary gluonic distribution (24) is much more 
concentrated at small x than that used in ref.[18] - xG 7 (x) ~ x 5 (l — x) a , 
where a = 0.1 at normalization point Ql = 0.3GeV 2 and a ~ 1 at Q 2 ~ 
5 lOGeF 2 ; iii) the difference in the boundary conditions for valence quarks 
- eq.(17). 

The change of the normalization point to Ql = 2GeV 2 (instead of Ql = 
lGeV 2 ) decrease gluonic distribution in our calculation by 10-12% (Fig. 8b). 
Therefore, the Ql dependence of the final gluon distribution is weaker in the 
variant II than in the variant I. 

The results of ref.[22] are shown by dashed line in Fig. 13, and results of 
[23] -by crosses. 

The negative gluon condensate contribution to the boundary condition for 
valence quarks (eq.17) results in flattering of G{x) at small x ~ 0.02 — 0.05. 
Result of fit, made in [20], are higher than ours (in all six choises of their 
input) at low x, because their boundary condition for G(x) is 2.5 5 times 
larger, than ours. In Fig. 14 we compare our results (II variant, H2) whith 
three sets of boundary conditions, offered by [20]. (Three other sets are even 
more larger and we don't show them). 

In the considered above two variants of the boundary conditions (17,18) 
and (18,24) the gluonic distributions in the real or weakly virtual photon 
(I p 2 |< O.bGeV 2 ) are essentially different. It is impossible to remove this 
difference by varying the normalization point Ql in the first variant of the 
boundary conditions. The experimental investigation of the gluon distribu- 
tion in photon can shed light on the problem if there are or not intrinsic 
gluons in the photon. 

The calculation of the whole F<i structure function can be done in a sim- 
ilar way. Like the gluon distribution, quark distribution functions are given 
by the sum of partial solutions of inhomogeneous equations with boundary 
conditions q mh (x, Qq) = plus general solution of homogeneous equation 
q hom (x,Q 2 ) with boundary conditions, following from (1). As we are in- 
terested only in the terms, proportional to In Q or a s ln Q, so the gluon 
contribution and a s -corrections to quarks contributions in (11a) should be 
calculated in In Q 2 order. This means, that this contributions may be ac- 
counted only in the solution of inhomogeneous equations since the terms, 
proportional to In Q 2 , arises only there. So, if we write ^(x) in the form 
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F 2 (x, Q 2 ) = F™ h (x, Q 2 ) + F^x, Q 2 ), 

where 

Ft\x, q 2 ) = x\ [ e e 2 (r"(^, q 2 ) + ^ ^ © c^, q 2 )+ 

+ <e 2 >^ 5 G 0G(x,Q 2 )| 

then 

F 2 ^(x,Q 2 ) = * £ e 2 ^(x,Q 2 ) 

i=l 

and the boundary conditions for q hom (x, Qq) are the same as (17) for the I 
variant (the II variant gives almost the same results and we will not discuss 

it). 

The results for F 2 are shown in Figs. 15-18. (The notations are the same 
as for the gluon distribution function, Figs. 4-7). The results for the real pho- 
ton case (thick solid line in Figs. 15, 16) are compared with the experimen- 
tal data of T PC /2-y [24,25], PLUTO [26-28] and OPAL [29] collaborations. 
For comparison the results of [18] (dashed line) and [20] (crossed line) are 
shown. One can see, that our results are in a good agreement with the data 
at 0.2 < x < 0.7 and Q 2 = 5 + 4lGeV 2 where our approach is valid. At 
x < 0.2 the gluon condensate contribution becomes too large and our results 
are unreliable. We should notice that, of course, at Q 2 > 40 -i- bOGeV 2 the 
role of c-quarks is not negligible at x < 0.3 and one has to take them into 
consideration. 

The results of the structure function F 2 calculations for virtual photon are 
presented in Figs. 17,18. The thick solid line corresponds to the transversally 
polarized photon case, the dashed one - to the unpolarized virtual photon 
case (the latter will be discussed below). 

Note that our results comparatively weekly depend on the renormalization 
scheme (MS or DIS 7 ), as it can be seen from Fig. 19, where our results are 
plotted in these two schemes. An approximate scheme independence arises 
because the scheme difference in the solution of inhomogeneous equation (see 
[16]) is compensted by the opposite difference of the solution of homogeneous 
equation. The main reason for this circumstance is that we have no separate 
boundry conditions for homogeneous equation (which is independent of the 
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solution of inhomogeneous equation and is determined from any physical 
ideas like VDM), but have boundary condition for total F 2 (x, Qq). 

Finally, let us discuss the nonpolarized virtual photon case. It is reduced to 
trivial redefinition of boundary conditions for homogeneous equation (17,18). 
Here we now should use F 2 = F 2 + (1/2) F L , F 2 being defined in (1) and 
F 2 can be found in [1] 

F 2 L = ~ Ee'^l-s) ^ 



7T " ' V ' n 2 !^L_ 

r x(l-x) 

m q - is the quark mass. 

Note that there is no contribution from gluon condensate in F 2 (see [1]). 

The results for the unpolarized virtual photon case are about 20% higher 
than for transversally polarized one and are shown by dashed line in Figs. 5-7, 
10-12 for gluon distribution and in Figs. 17,18 for F2 structure functions. For 
comparison, the results of [19] (disregarding any heavy quark contribution) 
are shown in Figs. 10, 11, 18 by crossed line. 
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Figure Captions 



Fig. 1. Bare Diagram. 

Fig. 2. Hadronic contribution into structure function. 

Fig. 3. Some diagrams corresponding to quark and gluon direct 
production. 

Fig. 4. Gluon distribution in the real photon (thick solid line); 

thin solid line and dotted line - respectively, the con- 
tributions from hadronic part and from gluonic conden- 
sate, the latter taken with opposite sign. The results are 
given for Q 2 = 5, 10, bOGeV 2 . The boundary conditions 
eq.s (17,18) - I variant, correspond to the normaliza- 
tion point Ql = lGeV 2 . The dashed line stands for the 
result of ref.[16]. 

Fig. 5. Gluon distribution for virtually polarized photon (p 2 = 
-0.5GeV 2 ,Ql = 1 GeV 2 )] thick and dashed lines, re- 
spectively - transversally polarized and unpolarized vir- 
tual photon case, hadronic part and module of gluon 
condensate contribution are shown by dashed and dot- 
ted lines, respectively. 

Fig. 6. The same as Fig. 5 for p 2 = -lGeV 2 ] Ql = 
4GeV 2 ; Q 2 = 50GeV 2 . 

Fig. 7. The same as Fig. 5 for p 2 = —2GbVq\ Ql = 
6GeV 2 ; Q 2 = WOGeV 2 . 

Fig. 8. Gluon distribution function in the real photon at Q 2 = 
lOGeV 2 and different Q 2 , Q 2 = lGeV 2 - thick solid line, 
Ql = 2 - thin solid line), a) Boundary conditions eqs. 
(17,18) - I variant, b) Boundary conditions eqs. (18,24) 
- II variant. 



23 



Fig. 9. Gluon distribution in the real photon for the case of the 
boundary condition (24) and (18) for quarks - II variant 
of the boundary conditions, Q\ = lGeV 2 . The notations 
are the same as in Fig. 4. 

Fig. 10 The same as Fig. 5 for the II variant (eqs. 18,24). The re- 
sults of [19] are shown by the crossed line for comparison. 

Fig. 11 The same as Fig. 6 for the II variant (eqs. 18,24). The re- 
sults of [19] are shown by the crossed line for comparison. 

Fig. 12 The same as Fig. 7 for the II variant (eqs. 18, 24). 

Fig. 13 Gluon distribution in the real photon, Q 2 = lOGeV 2 . 

Thick solid line-our results (II variant), thin solid line - 
our results (I variant), dotted line - the results of [18], 
dashed line - the results of [22] , crossed line - the results 
of [23]. 

Fig. 14 Gluon distribution in the real photon; thick solid line - 
our results, thin solid line, dotted and dashed lines - the 
results of fit [20]. 

Fig. 15 Structure function of real photon for three flavours 
(thick solid line) at Q\ = IGeV 2 for Q 2 = 5 -r IbGeV 2 . 
Thin solid line and dotted line contribution from 
hadron part and gluon condensate, the later - with op- 
posite sign. Experimental data (Refs. [24-29]) are shown 
for comparison. Also results of Ref. [18,20] are shown by 
dashed and crossed lines respectively. (On Fig. 15a re- 
sults of PLUTO collab. are at 5.3GeV 2 and TPC - at 
b.lGeV 2 .) 

Fig. 16 The same Fig.15, for Q 2 = 20, 45GeV 2 . 
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Fig. 17 Structure function of transversally polarized virtual pho- 
ton (thick solid line) and unpolarized virtual photon 
(dashed line). Thin solid line and dotted line - as in 
Fig.15. p 2 = -0.35GeV 2 , Q 2 = IGeV 2 . 

Fig. 18 The same as Fig. 17 for p 2 = -IGeV 2 , - 2GeV 2 and 
Ql = 2GeV 2 , 4GeV 2 respectively. Cross line - the re- 
sults of Ref.[19]. 



Fig. 19 Comparison of our results in MS (thin solid line) and 
DIS 7 (thick solid line) schemes for gluon distribution 
and structure function in real photon case. 
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